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Abstract 



The linear and nonlinear processes in ferromagnetic films at low temper- 
atures (T <C T c ) are studied in a microscopic theory. Both the long-range 
magnetic dipole-dipole and the Heisenberg exchange interactions to nearest 
and next-nearest neighbors are included. The results obtained for the lin- 
earized spin-wave spectrum are compared with previous macroscopic theo- 
ries. For ultrathin films (or for large wave vectors) the microscopic theory 
provides important corrections. The nonlinear dynamics of the spin waves 
are studied through a finite-temperature perturbation theory based on Feyn- 
man diagrams. We obtain explicit results for the energy shift and damping 
(or reciprocal lifetime) of the dipole-exchange spin waves due to all possi- 
ble three-magnon and four-magnon processes involving combinations of the 
surface and quantized bulk spin waves at low temperatures. To investigate 
different dipole interaction strengths (relative to the exchange) numerical re- 
sults are presented using parameters for Fe, EuO, and GdChj. 
75.70.Ak, 75.30.Ds, 76.50.+g 
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I. INTRODUCTION 



Ultrathin ferromagnetic films are of great interest due to their widespread technological 
applicability. 1 This is because these films can have properties, both static and dynamic, 
that differ significantly from bulk samples of the same material. For example, it is well 
known theoretically that, in these quasi- two-dimensional systems, the short-range exchange 
interactions alone are not sufficient to establish a ferro magnetically ordered ground state, 2 ' 3 
and it is necessary to take into account the anisotropy and long-range character of the 
dipolar interactions. Regarding the spin dynamics in ferromagnetic films, a macroscopic 
(or continuous medium) theory for the dipole-dominated regime was given by Damon and 
Eshbach 4 in terms of magnetostatic modes. They identified a surface branch to the spectrum, 
now known as the Damon-Eshbach (DE) mode. Subsequently there has been extensive work 
to generalize the magnetostatic theory to include also the exchange effects, leading to dipole- 
exchange theory (for reviews see, e.g., Refs. 5-7). In these theories the DE surface mode 
becomes modified and various quantization effects become evident for the bulk modes. In 
addition, Benson and Mills 8 considered a microscopic formalism for the dipole interactions 
in a thin ferromagnetic film, obtaining a rapidly converging series for the dipole sums. More 
recently, their approach has been extended to a description of the dipole-exchange spin waves 
(SW) in various ferromagnetic and antiferromagnetic films. 9-11 

All of the studies mentioned above neglected the higher-order effects that give rise to 
SW interactions. These occur because the SW are not exact eigenstates of the magnetic 
Hamiltonian and they are the source of several nonlinear effects, 6 ' 7 such as parallel pumping 
of SW, auto-oscillations and the transition to chaos, soliton formation, and energy renormal- 
ization and damping of the SW modes. Previous studies of the SW interactions in films were 
carried out from a macroscopic point of view, 5-7 where the dipolar terms were introduced 
using Maxwell's equations. Again, this should be appropriate for small wave vectors and for 
moderate film thickness. 

In infinite ferromagnets, either with or without the inclusion of the dipolar interactions, 
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the treatment of the interactions between the bulk SW modes is well established. 12 The dom- 
inant contributions at T <C T c come from three-magnon confluence and splitting processes 
and from four-magnon scattering processes. In the special case of a Heisenberg ferromagnet, 
where there are only exchange interactions, the three-magnon process is absent. However, 
for a finite ferromagnet, the presence of surfaces leads to a richer SW spectrum that may 
consist of one or more localized surface SW modes as well as the bulk modes (which become 
quantized in a film geometry). The occurrence of surface and modified bulk modes in a 
finite ferromagnet will give rise to more complicated schemes or SW interactions, since the 
scattering processes may now involve both types of modes. 

There has previously been a number of calculations for SW interactions in semi-infinite 
Heisenberg ferromagnets at low temperatures T <^ T c using a Hamiltonian approach. In 
particular, results for the damping and the energy renormalization of the surface SW were 
obtained in some special cases by Tarasenko and Kharitonov 13 and Mazur and Mills. 14 More 
detailed results, which were applicable to both surface and bulk (volume SWs and induced 
pinning through surface anisotropy and modified surface exchange, were derived by Kontos 
and Cottam 15 ' 16 using a diagrammatic perturbation method. 

In the above-mentioned studies for Heisenberg systems, the dipole-dipole interactions 
were ignored, leading to a description of SW interactions at T <C T c only in terms of the 
leading-order four-magnon processes. When the dipolar terms are included there are also 
three-magnon processes involving the splitting of SW into two SW modes and the conflu- 
ence of two SW modes into a single SW. For example, Rahman and Mills 17 calculated the 
renormalized energies for the bulk and Damon-Eshbach surface SW in the dipolar-exchange 
regime. For ferromagnetic films in this regime of long wavelengths (or small wave vectors), 
there have been extensive calculations of the nonlinear SW properties using a macroscopic 
(or continuous-medium) type of theory. 6 ' 7,17,18 A nonlinear microscopic theory, based on a 
Hamiltonian approach, allows us to study the SW interactions and the renormalization of 
the discrete SW modes. 

The aim of this paper is to provide a microscopic theory to calculate linear and nonlinear 
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properties of the localized surface SW and the quantized bulk SW in ferromagnetic films. 
This is done for all wave-vector values and for any film thickness. The theory includes the 
effects of the short-range exchange coupling and the long-range dipole-dipole interactions 
between spins. The latter interactions are evaluated in a similar manner to recent work 
for ferromagnetic and antiferromagnetic films in the linear SW regime. 9-11 We obtain the 
dependence of the SW energies (for each of the discrete branches) on the number of layers of 
the film, the variation of the exchange parameters (including nearest and next-nearest neigh- 
bors), the strength of the dipolar interaction, and the in-plane wave vector. We also develop 
a perturbation formalism to study the linear and leading-order nonlinear SW processes of in 
ferromagnetic thin films at low temperature T <T C . We use a Hamiltonian formalism to 
find formal expressions for the SW interaction terms involving products of three and four bo- 
son operators, representing the three- and four-magnon processes respectively. The method 
basically involves diagonalizing the non-interacting part of the Hamiltonian that is bilinear 
in the boson operators by making a transformation to a new set of boson operators, and 
then using it to define "unperturbed" Green functions. The perturbation terms, involving 
either three or four boson operators, are each transformed to the new boson operators and 
used to define interaction vertices within a Feynman diagram formalism. 

This paper is arranged as follows. Section II describes the Hamiltonian with exchange, 
Zeeman, and dipolar terms, as well as the assumed geometry of the ferromagnetic film. 
The method basically involves generalizing some previous calculations for nonlinear SW 
processes in semi-infinite Heisenberg ferromagnets 15 to include the dipolar interactions and 
the effect of finite film thickness. The transformations of the dipole sums that enter into our 
Hamiltonian were calculated according to the method described in Ref. 8. Results for the 
linear (noninteracting) SW spectrum of the film are derived in Sec. Ill, where we include 
numerical applications to different materials. In Sec. IV we analyze the higher order terms 
in the Hamiltonian that lead to SW interactions. At low temperatures T <C T c , it is shown 
that the dominant interaction processes involve three- and four- magnon scattering, by 
analogy with infinite ferromagnets. 12 . In Sec. V we introduce the Green function formalism, 
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associated with the diagonalized Hamiltonian of the system, and employ it to obtain 
general expressions for the energy shift and damping of the magnons in a ferromagnetic film 
in terms of a SW self-energy contribution. A diagrammatic perturbation technique is then 
employed in Sec. VI to calculate the SW self-energy due to the three- and four-magnon 
processes. In Sec. VII we discuss the numerical results applied to three different materials: 
Fe, where the ratio of dipolar to exchange strengths is very small; GdCl 3 , where the dipolar 
terms are relatively strong; and EuO, which represents an intermediate case. The overall 
conclusions are in Sec. VIII, while some of the mathematical expressions are quoted in the 
Appendices. 

II. HAMILTONIAN AND FILM GEOMETRY 

We consider a ferromagnetic film with N atomic layers of spins (with quantum number 
S) arranged on a simple cubic lattice with lattice constant a. The external magnetic field 
H Q and the static magnetization M are assumed to be parallel to z direction, the surface 
of the film is in the x — z plane, while the y-axis is perpendicular to the surface. We can 
generalize to the other crystal structures, surface orientations, and magnetization directions 
in a straightforward manner. 

The system will be represented by the Hamiltonian 

H = -\Y. JiA ■ - 9I*bH £ s; + \{g^ B ) 2 £ DgS?Sf (1) 

ij i ij 

where Sj is the spin at site i, and is the exchange coupling between sites labeled i 
and j. It is assumed that the exchange coupling is J\ and J2 for nearest and next-nearest 
neighbors respectively, and zero otherwise. The second term of Eq. (1) describes the Zeeman 
interaction of the spins with an externally applied field H . The dipole-dipole interaction 
between spins in different lattice sites is represented by the last term in Eq. (1), where a 
and P denote components x, y, or z, and 



where = — r*, and the case % = j is excluded from the sums in Eq. (1). Here we 
consider, for simplicity, that the exchange constants do not change near the surface of the 
film, and surface anisotropy fields are also neglected. However, we can easily include these 
changes in our calculations. 

The Hamiltonian may be written in terms of boson operators using the Holstein- 
Primakoff transformation, 12 and can be expanded as H = + + + ... (apart 
from a constant), where ff ( m ) denotes the term involving a product of m boson operators. 
The non-interacting (linear) SW modes are described by the quadratic term that has 
the following form 

if ( 2 ) = {A { nl> (qK„aqn' + #2 CO) [ a qna_ q n' + <4»aL qn '] } ■ (3) 

qnn' 

Here we are using a representation of the boson operators a* and a in terms of a 2D wave 
vector q = (q x , q z ) parallel to the film surfaces and indices n and n'{— 1, 2, . . . , N) that label 
the atomic layers parallel to the surface. The amplitude factors are 

A®,(q) = LiJLBH + S 



U n (0) -Un(q) +Vn,n+l(0) +U n , n _i(0) - (<^ B ) 2 £ D%„ (0) 



+ (^r^(q)L (4) 

= ^ B ) 2 [I«q) - ^(q) - 2^„,(q)] . (5) 

Here we have introduced 

u n (q) = 2Ji<ri(q) + 4J 2( r 2 (q) (6) 

and 

Un,n±i(q) = Ji + 2J 2 <ri(q) (7) 

as intra-layer and adjacent-layer Fourier transforms of the exchange interactions, with 
<7i(q) = [cos(g x a) + cos(g z a)] and cr 2 (q) = cos(g x a) cos(g 2 a). The quantities f)"f, are analo- 
gous Fourier transforms with respect to q of the dipole interactions defined in Eq. (2). The 
expressions for these terms are given in Appendix A. 
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In the next section we discuss the linearized SW spectrum, which is obtained using 
H( 2 \ before proceeding in Sec. IV to a description of the nonlinearities produced by the 
higher-order terms. 



III. LINEARIZED SW SPECTRUM 

In the linear approximation the SW spectrum is obtained by applying the standard 
equation of motion ihdA/dt = [A, H] (for any operator A) to the boson operators a qn and 
a qn in layer n, with the Hamiltonian H replaced by H^ 2 \ Taking % — 1 and assuming a 
time dependence for the modes like exp(— iut), we obtain the sets of equations 

"< n = E K 2 i(qKn + p£'(q) + ^(-q)]^^} , (8) 

n' 

- ua^ nl = £{[£©(-<*) + sS(q)]a qn + 4fi(-q)a_ q „,} . (9) 

n' 

There are 2N coupled equations altogether, where N is the number of layers of the film. 
Using the symmetry properties A^,(q) = A^,(— q) and B$ n (—q) = B^,(q), which arise 
from the properties of the dipole sums (see Appendix A) and the property Jjj = Jjj, the 
condition for there to be non-trivial solutions of Eqs. (8) and (9) can be expressed as 



det 



0. (10) 



A( 2 )(q) -cul^v 2B( 2 )(q) 
2B( 2 )(q) A( 2 )(q)+^I W 

Here A^ 2 ^(q) and B^ 2 ^(q) are N x N matrices with elements defined in Eqs. (4) and (5) 
respectively, and Ijv is the N x N unit matrix. The discrete SW frequencies in the film 
correspond to the N solutions for positive u of Eq. (10). They are degenerate in magnitude 
with the negative-frequency solutions. 

Numerical applications of the above theory have been made using parameters appropriate 
to thin films of the ferromagnets GdCl 3 (T c w 2.2 K) and EuO (T c w 69 K). These are 
chosen as representing different strengths of the dipole interactions relative to the exchange 
interactions. They correspond, in fact, to increasing exchange strengths, as can be inferred 



from the increasing T c values. In general, the dipolar interactions will have their most 
significant effect on the long- wavelength (small q) SW properties in each material, while the 
exchange effects dominate at shorter wavelengths. We note that thin-film samples of EuO 
have been employed in SW experiments, e.g., using Brillouin light scattering and microwave 
resonance techniques. 19 ' 20 Consequently, their magnetic parameters are fairly well known. 
In the case of GdCl3 the approximate parameters can be deduced from a study of the SW 
spectrum in bulk samples (e.g., see Ref. 21). The ratio of AttM (which characterizes the 
dipolar strength) to the bulk exchange field H Ex is approximately 1.5 for GdCl 3 and 0.063 
for EuO. 

It is well known that, in Heisenberg ferromagnetic films, exchange-dependent surface 
SW modes may occur 5 ' 7 that have properties quite different from the Damon-Eshbach (DE) 
surface mode 4 ' 5 of the magnetostatic theory. They occur, for example, if the exchange 
coupling between spins in the surface layers is modified or if effects of next-nearest neighbor 
exchange interactions are included. Here we study such exchange surface modes (as well as 
the DE type of surface modes) by including the next-neighbor exchange interactions J 2 in 
our numerical calculations for ultrathin films. 

We first show calculations for a 16-layer EuO film, taking ArcM = 2.4 T, H Ex = 38 T, 
and Hq = 0.36 T. Note that the bulk exchange field is given by Q[XbHex = 6S(J\ + 2 J 2 ) for 
the model simple-cubic structure. Here J\ and J 2 are effective exchange parameters because 
EuO normally has a fee structure. In Fig. 1 plots are shown for the frequency in GHz 
(converted using 7 = 28 GHz/T) of the lowest few discrete SW branches versus the ratio 
J2/J1 and versus q x a/n (which ranges from at the Brillouin zone center to 1 at the zone 
boundary) for small q x . We have taken the case of q z = (the Voigt geometry), for which it is 
known that there is a Damon-Eshbach (DE) surface mode 4 in the magnetostatic continuum 
limit where exchange is neglected. In Fig. 1(a) the dependence of the two lowest branches 
of the SW dispersion relation on J2/J1 is shown for two different values of q x a/n. For q x 
we note that the first branch has almost no dependence on J 2 / ' J\. That is expected since 
this lowest branch is essentially the DE mode. The second branch has a smooth dependence 
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on J 2 . For q x = 0.01 both modes increase their energy as J2/J1 increases. These effects 
are more noticeable in Fig. 1(b), where a plot of the frequency against the wave vector is 
shown for three different values of J-ij J\. The purely magnetostatic DE mode would tend to 
the approximate limiting value of 43 GHz with increasing q x . While it can be seen in both 
cases that there is a relatively flat (dispersionless) branch near this frequency, the upward 
curvature apparent at larger q x is due to the exchange terms. In the region near 35 GHz 
for small q x a/ir where two branches (the lowest bulk mode and the analog of the DE mode) 
come close together, but do not cross. This phenomenon of a near "crossover", with mode 
repulsion, signifies a strong mixing (or hybridization) of two modes. 

In Fig. 2 we show the behavior of the lowest four SW branches for very small wave 
vectors in a GdCl 3 film with iV = 16, considering J 2 /Ji = (dashed lines) and J 2 / ' J\ = 0.25 
(solid lines). This shows the effects of the next-nearest neighbor exchange. For this material, 
we use the parameters AirM = 0.82 T, H Ex = 0.54 T, and H = 0.36 T. It can be seen that 
a purely magnetostatic DE mode starts at q x = at the frequency 18.25 GHz in both cases. 
By contrast with the same film thickness for EuO in Fig. 1, we can see that the DE mode 
in the GdCl 3 case affects several of the bulk modes of the spectrum, which is a consequence 
of the larger AirM/H Ex in this material compared with EuO. 

In the numerical examples shown above, the limiting frequency of the DE mode (as q x — > 
0) is very close to the value [H (H + AttM)} 1 / 2 predicted by the macroscopic magnetostatic 
theory 4 in the Voigt geometry. On the other hand, the magnetostatic theory predicts that the 
DE branch becomes flat, tending to a limiting value of Hq + 2ttM when q x is large compared 
with the reciprocal of the film thickness. This behavior is followed approximately in our 
theory, but there is an additional upward curvature to all SW branches due to exchange. 
This is in addition to the hybridization effects and the microscopic corrections to continuum 
dipole-exchange theories. 
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IV. HIGHER ORDERS: SW INTERACTIONS 



The higher order terms in the expansion of the Hamiltonian, i.e., the terms with 
m > 3, correspond to the nonlinear aspects of the SW dynamics including interactions 
between the SW modes. As in the case of bulk (effectively infinite) ferromagnets, 12 the 
leading-order effects for T <C T c come from the terms that have three and four boson 
operators (# (3) and H^). 

The term depends only on the dipolar interactions and has the form: 

^ = \ H A nl>( k ) [4n'«q-kn a qn + a q „a q -kn<w] , (11) 
kqnn' 

with the amplitude term given by 

(3) fU\ _ . /oC„2.,2 mxz 



4»(k) = V2 V/4 TO(k) - iDl z n ,{\,)} . (12) 

This represents the splitting and confluence processes inherent in the three-magnon pro- 
cesses. 

By contrast, the four-magnon interaction has parts involving both dipolar and 

exchange terms: 

H(4) = 2 E {^itl'( k )4'n4-k'n a q-knOkn' +^itl'(- k >L'4-k'n a q-kna k n 
kk'qrm' 

~~ 2-B^'(k — k^a^/a^k/jjOq-knakn' + C'in'( k )[ a qn a< i- k - k ' nak ' riakn ' + h.c\} . (13) 



In this case the amplitude coefficients are defined by 

1 



4 4 i(q) = \ 



+ fn,n+l(q)^n',n+l + ^n',n-l(q)£n',ra+l + ^/i B ) 2 Z^ n ,(q) 



(14) 



a$(q) = 4 4 i(q) - ^) 2 ^(q), (is) 



cffi(q) = -(i/5)BS(q). (16) 

This interaction represents the scattering of a pair of magnons into another pair of magnons. 
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We note that each term in Eqs. (11) and (13) conserve the 2D in-plane wave vector, since 
there is translational symmetry in the x — z plane. The dependence in the y-direction is 
taken into account the the summations over layer indices n and n' . The above H^> and 
interactions are the analogs (in a microscopic theory) of the expressions used in macroscopic 
theories to describe a wide range on nonlinear effects, 5,6 as mentioned earlier. For example, 
they describe processes whereby energy from a given SW (with, say, wave vector q and from 
the discrete branch v) can be transferred to other SW modes of the film. Thus they can 
be used to calculate the damping (or reciprocal lifetime) of any given SW, as well as the 
shift (or renormalization) in its energy. This can be accomplished using a diagrammatic 
perturbation method with Green functions. 



We establish a Green function formalism for the film geometry with iV atomic lay- 
ers parallel to the surfaces. The method is analogous to that used in previous diagram- 
matic perturbation calculations for the effects of SW interactions in semi-infinite Heisenberg 
ferromagnets, 15,16 but generalized here to include the dipolar terms and the quantization of 
the modes due to finite film thickness. 

The starting point to the diagrammatic formalism is the expansion of the Hamiltonian 
into terms containing products of two, three, and four boson operators. The relevant ex- 
pressions for these terms are given in Eqs. (3), (11), and (13). The bilinear term H^, which 
we shall treat as the non-interacting part of the Hamiltonian, can first be rewritten as 



V. THE GREEN FUNCTION FORMALISM 



= -- 



Tr [A (2) (q)]+-E^U(q)^ 



(17) 




where 



/ A 



A( 2 ) (q) 2C( 2 ) (q) 
2B( 2 ) (-q) A( 2 ) (-q) 



\ 



*(q) 



(18) 



V 



/ 



Here we have defined operators .4.* and *4 q that create and destroy wave vector q by 
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At 



y a_ q/ 



^4q 



\ a -q/ 



(19) 



respectively, where a q and a q are iV-component column matrices whose nth elements are 
a qn and a qi „. The tilde denotes a matrix transpose, and A^(q) and B^(q) are the N x N 
matrices with elements given by A^,(q) and B^,(q) defined in Eqs. (4) and (5). 

The first term in Eq. (17) is just a constant, while the second term provides the non- 
interacting (linear) SW spectrum. In fact, the eigenvalues of x(q) are the discrete SW 
energies E^ v where v = 1,2, ...,N is a branch label (each eigenvalue will occur twice 
in x(q))- Results for E^ u have already been obtained using the equivalent determinantal 
condition in Eq. (10). We now diagonalize by making a linear transformation to a new 
set of boson operators « q and a q , which satisfy the usual commutation relations and are 
defined by 



fl q,™ = J2 S nl (q) «q,« + T nl (q) a -q,H 
I 

a ln = S nl (q) + T nl (q) "-q,!- 



(20) 
(21) 



The S n i and T ni coefficients can be found by a straightforward generalization of a procedure 
due to White et al, 22 and the result is as follows. Suppose we rewrite the above transforma- 
tion in matrix form as A^ = 5**A q and A^ = SqA q , where A q and A q are defined similarly 
to A^ and Aq but in terms of the new operators. The jth column (denoted by Sj tq ) of the 
2N x 2N transformation matrix 5 q is then found by solving for the eigenvector in 

( A( 2 )(q) 2C( 2 )(q) \ 

. Si,q = ±Ci (q) S,,q, (22) 

^ -2B< 2 ) (-q) -A( 2 ) (-q) ) 

where the + sign is taken for j e 1, . . . , N and the — sign for j e N + 1, . . . , 2N. The H^> 
part of the Hamiltonian then assumes the simple form (apart from a constant term): 



ff (2) = E^( 1 + 2a W)- 



(23) 
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To establish the diagram technique we now introduce aiVxJV matrix G(q,iu m ) with 
matrix elements defined as the causal Green functions involving the transformed boson 
operators: 

G w /(q,iw m ) = ((a q „;aL)). . (24) 

in a conventional notation. 13 Here iu) m = 2nmi/j3 is an imaginary boson frequency, where 
(3 = 1/ksT and m takes all integer values from — oo to +00. When these Green functions 
are evaluated using just in Eq. (23) for the Hamiltonian the resulting non-interacting 
Green functions are diagonal in the labels v and v' and are given by 

G° uu ,{q,iuj m ) = Gl(q,iu m )5 u y = (~) - ^— S »y- (25) 

On making an analytic continuation iu m — > uj + i0 + to real frequency u, it is seen 
that there is a simple pole at the values E^ v . These correspond to the non- interacting 
excitations, comprising the discrete bulk and surface SW modes, which we have already 
been in the linear case. In the higher orders of perturbation, i.e. when and 
are included, the "interacting" (or renormalized) Green functions will have modified poles, 
corresponding to a renormalization in energy by an amount AEq jV and a damping r q;/ . 
Also, in general, the interacting Green functions will have off-diagonal terms in the v and 
v' labels. As in standard diagrammatic formulations for interacting boson systems, the 
interacting and non-interacting Green functions are related to one another by means of 
the Dyson development operator. 23 This quantity can be expanded as an infinite series in 
powers of the higher-order Hamiltonian parts. The terms in the expansion can conveniently 
be evaluated using a diagrammatic representation with proper self-energy terms H vv i (q, iu m ) 
introduced to account for the interactions. The choice of the self-energy terms within the 
perturbation expansion depends on the explicit form of and and will be discussed 
in the following sections. 

First, the non-interacting Green function (^^(q, iou m ) will be drawn as a solid directed 
line as shown in Fig. 3(a). There are then four possible types of proper self-energy terms, 
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as defined in Fig. 3(b), depending on whether the G° lines are entering or leaving. All of 
these may be nonzero when there are magnetic dipole-dipole contributions to the Hamilto- 
nian. The renormalized Green functions are given by a generalized Dyson series, which is 
represented diagrammatically in Fig. 4. This takes the algebraic form 

Gw(q,iu; m ) = G°(q,^>m)<W + G°(q, iu; m )£+7(q, «w m )G°,(q, iu m ) 

+ H G °(q ; ^m)£w",(q, ^m)G°„(q, iu; m )E+„~'(q, ^ m )G°,(q, iu; m ) 
q, ict> m )E^y) 

We note that the only term linear in the self-energy parts is the one involving E + ~. The 
terms with E ++ and E occur in diagrams quadratic (and higher) in the self-energy parts, 
and E~ + occurs only in diagrams that are cubic and higher. If we ignore, for the moment, 
the effects of E ++ , E , and E _+ , the infinite Dyson equation series for Gw(q, iuj m ) is 
summable, and the formal solution can be written in a matrix form as 



G(q, iu m ) = G°(q, iu m ) + G°(q, iu m )'S + (q, iw m )G(q, iw m ). 



(27) 



The solution of this matrix equation is 



G(q,iw m ) = \ fG°(q,iw m )l - E + (q,iw. 



-i 



(28) 



By solving for the complex poles of this interacting Green function we may deduce the 
renormalized SW energies AE^ V and their damping r qi ,. The condition is 



det 



G°(q, iu m )] - E + (q, iu m ) \ = 0. 



(29) 



After substituting for the matrix elements of [G°(q, iu m )] 1 using Eq. (25) this becomes 



# q ,i - ium - ^fi (q, iu m ) ^T,f 2 (q, iuj m ) 

^Ejf (q, iw m ) £ q ,2 - lW m - JS22~(q> ««m) 



0. 



(30) 



At low temperatures T <C 7^, where the perturbative expansion is expected to be good, 
an approximation to Eq. (30) (and hence to AE and T for each SW branch) may be obtained 
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by expanding the determinant to lowest order in the off-diagonal elements. The results may 
eventually be expressed as 

A£ q ,„ + iT^ = -^E+-(q, £ q ,„ + *0+) + 0[E+7(q, £ q ,„ + *0+)] 2 , (31) 

where v' ^ v. In deriving Eq. (31) we have assumed that |AE + ir| <C E, and so the matrix 
elements are evaluated there at the unrenormalized SW energy. 

When the other types of self-energy terms (i.e., S ++ , £ , and E _+ ) are taken into 
account, it is found that they modify Eq. (31) only by adding terms of second order (and 
higher) on the right-hand side. Therefore, to leading order (neglecting the quadratic self- 
energy terms) we have 

A£ q ,„ ~ ~Re E+-(q, E^ + t0 + ), (32) 

r q ,„ ^ -^Im E+-(q,^ + 20 + ), (33) 

for the energy shift and damping, respectively, of SW branch v at in-plane wave vector q. 
In the following sections the above expressions are used to study the effects of the three- 
and four-magnon interaction processes. 

VI. THREE- AND FOUR-MAGNON PROCESSES 
A. Three-Magnon Processes 

The term in the expansion of the Hamiltonian is given in Eq. (11). It arises entirely 
from the dipole-dipole interactions and is of third order in the boson operators. It describes 
splitting processes, in which one magnon is absorbed and two are emitted, and confluence 
processes, in which two magnons are absorbed and one magnon is emitted. These can lead to 
transfers of energy between the different SW branches (surface and quantized bulk modes). 

The first step is to rewrite in terms of the new boson operators a q and a q using 
Eqs. (20) and (21), so that the diagrammatic formalism may be employed. The result is 
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H(3) = \ H a k,/i a q-k,Z 2 «q,i3 + V 2 ai )h Q! q -k,i 2 Q!k,i3 

/lZ2^3kq 

+ V 3 "M^q-k^^-q^ + V * «-q,Zi«q-k,Z 2 «k,« 3 , (34) 



where the amplitude factors Vi (with i = 1,2,3,4) are listed in Appendix B. The leading- 
order diagrammatic contributions to the proper self-energy term as required for the 
renormalization of any SW branch I, are shown in Fig. 5. Each of the dotted lines represents 
a V{ (as an "interaction vertex" ) according to the direction of arrowing on the three Green- 
function lines entering or leaving it (e.g., a V 4 vertex has three lines leaving). All of the 
diagrams in Fig. 5 are of second order in the vertices. Using the formal rules of diagrammatic 
evaluation for interacting boson systems, 23 we obtain the self-energy expression as 

n°(E qV 0+n o (^ q - q ^) + l 



^(0,^) = -^ {(Wa + Wt fY*'^ +W C - 
q'i'I" I tj ^ 1 " 



| w n\E^ v )+n\E^,,) + l | w f{Eq',y) - n (E q+q ,, t ,Q 1 (35) 

Eq-q',l" + -EqV — -Eq+q',Z" — Eqiji — lUJ m J 

The weighting factors W a ^,..., e are given in Appendix B and n°(E) = [exp(@E) — 1] _1 denotes 
the Bose-Einstein thermal factor at energy E. Next, using Eqs. (32) and (33), together with 
the well-known relation 

p(-)t^S(x), (36) 



(x ± iE) \x 

where P indicates that the principal value is taken in a summation over x, we obtain 

n°(£ q 'j') + n°(£_ q _ q ',,») + 1 



= peK + ^ 6 )^zl + Wc 

ci'l'l" { ^iJ" 



E-q-q'J" + Eqiji + -E qj / 



, w n°(E^)+n°{E^,) + l | T ^ n (g q >, y ) - n°(g qW ) ) (3?) 

Eq-q'J" + Eqiji — Eqi Eq+qi^i — -Eq',Z' — £q,Z J 

for the energy shift, and 

q'Z'Z" 

+ W d [n°(£? qf) ,0 - n°(£ q+q ,^)] <K# q + q ',z" - £ qV , - £ q ,z)} , (38) 

for the damping. The two delta functions in Eq. (38), which conserve energy and in-plane 
wave vector, correspond to the three-magnon splitting and confluence processes respectively. 
The applications of the above results for AE and T are discussed in Sec. VII. 
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B. Four-Magnon Processes 



Following the previous section we first rewrite the interaction Hamiltonian, in this case 
the four-magnon contribution given by Eq. (13), in terms of the new boson operators. 
After a straightforward but lengthy calculation we obtain the following expression: 

h(4) = 2 51 { A i a k'z 1 a q-kv 2 «q-k,* 3 «k,u + A 2 a qh o; q -k-k',i 2 Q;k',i3«k,k 

ilWsUkk'q 

+A 3 aj c ' i a k ',/ 2 Q; q-k-k',/3 Q; q,U + A ^ a lh a l-k',i 2 a ^-ih a -Kh 
+A 5 a_k'iia k '-q,/ 2 aq- k ,/3 a k,/ 4 } , (39) 

where the amplitude factors Aj {i = 1,2,3,4) associated with each operator term are given 
in Appendix C. These terms define the interaction vertices of the diagrammatic formalism. 
By contrast to the three-magnon case, the four-magnon interactions involve dipolar and ex- 
change terms and contribute even in the Heisenberg limit (see, e.g., Kontos and Cottam 15 ' 16 
for the case of a semi-infinite Heisenberg ferromagnet). 

According to Eqs. (32) and (33), the formal results for the SW energy shift and damping 
related to arise from the real and imaginary parts, respectively, of an appropriately- 
chosen self energy term in the Aj vertices. The leading-order contribution (denoted by E^) 
to the proper self-energy S + "(q, iu m ) is of first order in Aj and corresponds to the diagram 
in Fig. 6a, together with topologically similar diagrams. Again the dotted line represents an 
interaction vertex corresponding to the number of lines entering or leaving. On evaluating 
the diagrams we obtain the contribution to the energy shift of SW branch I as 

A£ qii =^e a n°(%) (40) 

q'Z' 

where the weighting term a arises from a combination of Ai vertices and is defined in 
Appendix C. It is the analog, for the thin film with exchange and dipolar coupling, of the 
Dyson SW interaction term for the case of bulk SW in an infinite Heisenberg ferromagnet 
(see, e.g., Ref. 12). 

It turns out that the self-energy part is real, and so it gives a vanishing contribution 
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to the damping. Therefore, to calculate the damping, it is necessary to consider the self- 
energy contributions Y,^ that are of second order in the Aj vertices. This corresponds to 
diagrams with a more complicated loop structure (and hence with more internal variables). 
In Fig. 6b we show just those diagrams for £( 2 ) that give rise to a damping term (i.e., those 
that have an imaginary part after the analytic continuation iu m — > Eqj + i0 + is made. The 
final result for the damping of SW branch / is 

IV = -T E { Q b K(#q'-q,i') K(£q'V") + l] " n°(£ q ,_ q «,,»') 
q'q"l'l"l'" 

X [^"(E'q//^") — 7l°(.Eq/_ qj j/) ]} 5 (£q",J" + -Eq'-q"^'" — Eqi-qji — Eq^) 

+6 C {n°(f? q « ,«)n°(S q ,_ q « )JW ) + [n°(f? q «,,«) + 1 

X [n^-Eq",*") + ^(Eqi-qii^/ii) + lj | 5 (Eqiijh + Eq-qiji + Eqt-qiijm — £q,z) 

+e d {n°(£q,_q, r ) x [n°(£ qV + l] - n°(f? q ,_ qi ,0 

X [n°(_E q /y//) - 7l (£q/_ qjJ /) ]} 5 (-Eq'+q",/'" — £q",Z» — Eq'-q,J' ~ ^q,*)} ■ ( 4 1) 

The second-order weighting factors Ob, c ,d, which multiply the three different delta functions 
conserving energy and in-plane wave vector, are given in Appendix C. 



VII. NUMERICAL RESULTS AND DISCUSSION 

Numerical calculations for particular materials are now presented based on the formal 
expressions for AEq^ and r qj j in the previous section. This provides us with predictions of 
the overall dependence on q and T for the different SW branches and allows us to compare 
the relative importance of the three- and four-magnon processes in various situations. 

We start by considering ultrathin films of EuO, where the dipolar interactions have 
a significant effect on the long-wavelength linear SW properties and the exchange effects 
dominate at shorter wavelengths (see Sec. III). We assume the same values for AttM and 
Hex as before, but for simplicity we ignore the next-nearest exchange J 2 . In Fig. 7a we show 
the damping of the lowest SW branch (labeled 1) for N = 8 and for temperature T <C Te- 
as a function of the in-plane wave- vector component q x . The mode under consideration is 
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the analog of the DE surface mode at small q x which is then modified by the exchange at 
larger q x . In this case it is found that the dominant damping mechanism comes from the 
three-magnon term in second order, i.e., from Eq. (38). The delta functions in Eq. (38) 
give Eq i = Eqppqijn ± -E q ',/', in accordance with conservation of energy and in-plane wave 
vector. For a SW with in-plane wave vector q and branch /, these correspond physically 
to the splitting (upper signs) and confluence (lower signs) processes. Here the labels I' and 
I" may be the same as, or different from, /. Along with the total three-magnon damping 
contribution (solid curve) we show in Fig. 7a the contributions from some of individual 
(/', I") processes. It is seen that several of the inter-branch processes (i.e., where at least one 
of I' and I" differs from I) play an important role. In Fig. 7b we make some comparisons 
between the total damping for the SW branch 1 with N = 8 (solid curve) and N = 16 
(dotted line). In this particular case the damping has a similar behavior for both values of 
N, although the individual inter-branch contributions (not shown) are different. In the same 
plot we show the total damping for branch 2 in the case of N = 8 (dashed curve). This is 
slightly larger than the damping for branch 1, mainly because there are SW branches above 
and below it, giving more possibilities to satisfy the required conservation conditions. 

It is interesting to compare the damping results predicted for EuO with those for other 
materials. We consider Fe where the ratio of dipolar to exchange strength is smaller than 
in EuO, as well as GdCl 3 where the ratio is larger. For the case of a GdCl 3 film with N = 8 
layers, the total damping at T <C T c is shown in Fig. 8, where we use the same AttM and 
Hex as in Sec. III. This presents a rather different behavior, particularly for small q x , when 
compared to a EuO film (see Fig. 7a) and also to a Fe film (see Fig. 9). The parameters 
used for the latter case are AttM = 2.14 T and H Ex = 2140 T. The qualitative differences for 
the dominant three-magnon damping in the case of GdCl3 is mainly because the influence of 
the dipolar terms (including the role of the DE surface mode) extends over a greater range 
of wave vectors than is the case for EuO or Fe. This makes it easier to satisfy the energy 
and wave-vector conservation conditions for the damping at smaller wave vectors. 

As mentioned earlier, the three-magnon damping has contributions due to splitting and 
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confluence processes. The temperature dependence of these processes arise from the combi- 
nations of Bose factors associated with each delta-function term. We note, in particular, that 
the second term (the confluence term) on the right-hand side of Eq. (38) vanishes as T — > 0. 
This is as expected since there needs to be a thermally-excited SW to participate in the 
confluence. By contrast, the first term in Eq. (38) does not vanish in the zero-temperature 
limit, since it is always possible (provided the conservation conditions are satisfied) for a 
SW to split into two modes. The numerical calculations presented above were all for the 
low-temperature limit where the three-magnon splitting is dominant. In Fig. 10 we show 
calculations for the temperature dependence of the three-magnon damping, taking the case 
of a EuO film (with the same parameters as in Fig. 7a) and a fixed value of the wave- vector 
component (q x /n = 0.5). It can be seen that, at low enough temperatures, the splitting pro- 
cess (solid curve) dominates and increases slowly with temperature, whereas the confluence 
process (dashed curve) increases very rapidly with temperature. 

The four-magnon damping, which is given by Eq. (41), contains three different processes 
proportional to the weighting factors 0&, 6 C and 9^, together with their Bose factors. The 
delta functions provide for conservation of energy and in-plane wave vectors. The only 
non-vanishing contribution as T — > is the term proportional to 6 C , which describes a 
splitting of an incoming SW into three modes. It is proportional to the dipole-dipole sums 
(i.e., it vanishes in the exchange limit) and is generally less important than the splitting 
process in Eq. (38) discussed above. The other two terms in Eq. (41) contribute to the 
four-magnon damping when T ^ 0. The first term (proportional to O^) is the analog of the 
scattering term in the damping for Heisenberg systems. 12,15,16 The incoming SW scatters off 
a thermally-excited SW into two other SW modes. The term proportional to 6^ corresponds 
to the incoming SW scattering off a pair of SW modes into another SW mode. It arises 
specifically as a consequence of the dipolar terms in the Hamiltonian. 

We now turn to a discussion of results for the SW energy shifts, which are given by Eqs. 
(37) and (40) for the three- and four-magnon cases respectively. In a low-temperature limit 
T <^ Tc the dominant contributions come from the three-magnon processes, specifically 
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from the terms proportional to W c and Wd in Eq. (37). Numerical results for AEqj versus 
q x for EuO films in the same three cases as Fig. 7b are shown in Fig. 11. The main difference 
(compared to the damping) is that the inter-branch terms are much more important. This 
arises basically as a consequence of AE^i and r q ^ coming from the real and imaginary 
parts, respectively, of self-energy terms in the diagram formalism, and these may have quite 
different dependencies. For example, the imaginary parts have delta functions, which lead to 
peaks and structural features in the damping, whereas the summations in the real parts tend 
to give a smoother behavior for the SW energy shifts. Hence AE has a weaker dependence 
on q x in Fig. 11 than was the case for the damping. The larger AE value and its in- 
dependence for the thicker film (where the linear SW modes are closer together) is due to 
the inter-branch terms in the summations. 

At low but finite temperatures with T C T c the four-magnon contribution to AE^i be- 
gins to play a role, particularly for the lower-energy SW branches, because the Bose factor 
in Eq. (40) becomes nonzero. Some numerical calculations are shown in Figs. 12 and 13 for 
EuO films with N = 8. In Fig. 12 the four-magnon contribution to the energy shift for the 
lowest SW branch (I — 1) is plotted as a function of wave- vector component q x at a fixed 
temperature (T = 0.04 K). It is seen that the qualitative behavior is quite different from that 
for the three-magnon contribution (see Fig. 11). We note also that the (/^-dependence in Fig. 
12 is different from that of the four-magnon energy shift for a Heisenberg ferromagnet, 12 ' 15 ' 16 
which is proportional to g 2 . at small wave vectors. This difference, compared to Heisenberg 
systems, is partly due to the the dipolar terms in the factor a in Eq. (40) and partly due 
to the dipolar contribution to the SW energy gap. In Fig. 13 we compare the temperature 
dependence of the three-magnon SW energy shift (solid curve) with that due to four-magnon 
processes (dashed curve). The latter quantity increases rapidly with temperature and even- 
tually dominates. This is again for a EuO film, but for a fixed wave vector corresponding 
to q x = 0.5n. 



21 



VIII. CONCLUSIONS 



We have developed a microscopic (Hamiltonian-based) theory for the dipole-exchange 
SW modes and the interactions between them in ultrathin ferromagnetic films. We first 
used the theory to study the linear dispersion relation of the discrete SW in EuO and 
GdCl 3 , including the effects of next-nearest neighbor exchange. The results obtained with 
the microscopic theory are in good agreement with those of the usual macroscopic theories 
in limiting cases where there are a large number of layers and the wave vectors are very 
small. Outside of these regimes it becomes necessary to use a microscopic approach, as in 
this work. 

Secondly the theory enabled us to develop the leading-order three- and four-magnon 
interaction terms in the Hamiltonian for the film. These three- and four-magnon processes 
can involve modes from the different discrete SW branches (surface or bulk) in all combi- 
nations. Formal expressions were found for the energy shift and damping of each branch. 
Applications were made to EuO, Fe, and GdCl 3 . In the case where the three- magnon pro- 
cesses dominate, it has been shown that there are significantly different types of behavior 
predicted for GdCl3 compared to EuO or Fe. This happens because the influence of the 
dipolar terms extends over a greater range of wave vectors in the GdCl 3 . As a consequence, 
the energy and wave-vector conservation conditions for the damping are more readily sat- 
isfied at smaller wave vectors. Our method enabled us to carry out calculations also at 
higher temperatures, where the four-magnon processes play a more significant role. The 
above-mentioned influence of the dipolar interactions also makes an important contribution 
to these processes. The major effect, compared to the four-magnon energy shift in a Heisen- 
berg ferromagnet without dipolar interactions, can be observed in the q x dependence of the 
four-magnon energy shift. 

In order to probe the effects discussed here, inelastic light scattering would be an appro- 
priate experimental technique. 19 An interesting extension of our work would be to apply the 
nonlinear theory in an examination of SW instabilities in ultrathin films under conditions 
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of parallel or perpendicular "pumping" by a microwave field. 
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APPENDIX A: DIPOLE SUMS 

Here we list the expressions and properties of the dipole sums required for this work. 
They are calculated as in Ref. 8. 

The sums of dipole interactions between spins located in different atomic layers (i.e., 
with y = || (n — n')|| 7^ 0) have the following expressions: 

(&,?,) =4* £ inl + qx/2)2 e W (-2\y\ llm ), (Al) 

l,m=— oo llm 



oo 



D y n v n , (q x , q z ) = -4tt £ llm exp (-2 \y\ llm ) , (A2) 

Lm=— oo 



DZ n ,(q x ,q,) = 4* £ { — + ^ exp (-2 | y| llm ) , (A3) 

i,m=— 00 Tim 



(fc, &) = -i47rs 5 n(y) £ (ttZ + g x /2) exp (-2 |y| 7im ) , (A4) 

Z,m=— 00 



£>ni' (?*» ?*) = -i4irsgn(y) J2 i^ 171 + ^/ 2 ) ex P (~ 2 7«m) , (A5) 



(*,,,) = 4. £ (T ' + fe/2)(Tm + g ' /2) exp(-2| j/i7 ,„). (A6) 

l,m=— oo Tim 

In the expressions above we denote 7; m = [(7rZ + q x /2) 2 + (jm + ^/2) 2 ], and Ki(x) is a 
modified Bessel function of integer order % for any variable x. 
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For the dipolar sums involving spins in the same layer (y = 0) we have the following 
expressions for the diagonal terms (involving the same superscript): 

Dn* n (Qx: q x ) = ~2S X (q x , q y ) + S z (q x , q y ) , (A7) 
D y n V n (5x, Qz) = Sx {q x , Qy) + S z (q x , q y ) , (A8) 
D n,n (Qx, Qz) = ~2S Z (q x , q y ) + S x (q x , q y ) , (A9) 



where 



j^g oo oo 

Sx (Qx, Qy) = — J2 E (7vm + q z /2) 2 cos(q x x)[K 2 (2x\irm + q z /2\)}, (A10) 

" x=l m=—oo 



oo oo 

\2 



16 

Sz(qx,q y ) = —J2 E (nm + q x /2) 2 cos(q z z)[K 2 (2z\irm + q x /2\)}. (All) 

" 2 = 1 TYl= — OO 



Also for the off-diagonal terms: 

oo oo 

£>n?n(9*,&) = 16£ E (7rm + g z /2)|7rm + g z /2| sin (q x x)[K 1 (2x\irm + q z /2\)} , (A12) 



x=l m=— oo 

while £>^ n (g x , g 2 ) = D y n z n (q x , q z ) = by symmetry Finally, when q x = q z = 0, we have 



*C (o, o) 



47T 



2 



9 



1 + 24 jr m 2 K 2 (2a;7rm) 



m=l 



(A13) 



APPENDIX B: THREE-MAGNON AMPLITUDE FACTORS 



The expressions for the amplitude factor Vi (i = 1, 2, 3, 4) in Eq. (34), and hence for the 
vertices in the diagrammatic representation, are given by 

^(k, q\k, l 2 , l 3 ) = ]T {^(k)^, Jk)^ 2 (q - k)^ 3 (q) 

nn' 

+4 3 i(k)^ 1 (k)T4(k-q)T„ i3 (q) 
+4 3 i(-k)T:,J-k)^ 2 (q - k)^ 3 (q) 
+4 3 i(-k)T:,J-k)T4(k - q)T n , 3 (-q) 
+4 3 i(-q)^ 1 (-k)^ 2 (q - k)T n , Ia (-q) 
+4 3 i(q)^(k)T4(k - q)^ 3 (q)} , (Bl) 
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V 2 (k, q \h, Z 2 , l 3 ) = £ {4 3 i(k)^aq)^ 2 (q - k)5 B ^(k) 

tin' 

+4 3 i(k)T4(-q)T n , 2 (k - q)5 n , 3 (k) 
+4 3 i(-k)^ i (q)^ 2 (q - k)T n ,, 3 (-k) 
+A^,(-k)T: h (-c L )T nl2 (k - q)T„, 3 (-k) 
+AS(-q)T* ii (-q)^ 2 (q - k)T„, 3 (-k) 
+4 3 i(q)^Jq)^ 2 (k - q)^ 3 (k)} , (B2) 

V 3 (k, q \h,l 2 ,l 3 ) = Ei^'W^CkJ^Cq-kjT^Cq) 

nn' 

+4 3 i(-k)T n *, 1 (-k)T4(k - q)5^(-q)} , (B3) 
V 4 (k, qj/x, Z 2 , Z 3 ) = £ {4 3 i(k)T nil (q)^ 2 (q - k)S n ,, 3 (k) 

nn' 

+4 3 i(-k)^ 1 (-q)T„ i2 (k - q)T„, 3 (-k)} . (B4) 

In terms of the above quantities in Eqs. (35), (37) and (38), we can now define the 
weighting factors Wi (i = a, b, c, d, e) as 

W a = [V, (0, q']Z", l', l') + V l (q\ q'|Z', Z", I')] [V 2 (0, q| J, /, I") + V 2 (q, q|Z, I", /)] , (B5) 

Wb = [V 2 (0, q'|Z", I', V) + V 2 (q', q'|Z', Z", J')] (0, q| J, I, I") + V x (q, q|Z, Z", J)] , (B6) 

W d = [V, (q', q|Z', I", Z) + V 1 (q- q', q|Z", Z', /)] 

X [y 2 (q', q|Z', Z", Z) + V 2 (q - q', q|Z", Z', Z)] , (B7) 

W e = [14 (q', q + q'|Z', Z, Z") + V l (q, q + q'|Z, I', I")) 

x [F 2 (q', q + q'| J", Z, J') + V 2 (q, q + q'|Z", V, I)] , (B8) 

and finally W c = w±w 2 + w 2 W c2 + w{W c3 + 21^4, where 

= V 3 (q', -q|Z', Z", Z) + \/ 3 (-q - q', -q|Z", Z', Z) , (B9) 
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W2 = V 4 (q', -q|Z, I", I') + V 4 (-q - q', -q|/, /', I") , 



(BIO) 



W c2 = V 3 (q, -q'| J, I", I') + V 3 (-q - q', -q'|Z", /, Z') 

+V 3 (q, q + q'|Z, Z', Z") + V3 (q', q + q'|Z', Z, Z") , (Bll) 

^ c3 = V 4 (q, -q'|Z', Z", Z) + V 4 (-q - q', -q'|Z', Z, Z") 

+V 4 (q, q + qT, J', J) + V 4 (q', q + q'|Z", Z, Z') , (B12) 

W c4 = [F 3 (q, -q|/', I", I) + V 3 (-q - q', -q|Z", I', I) + 
V 3 (q', q + q|Z', Z, Z") + V 3 (q, q + q'|Z, I', I")} 
x [T/ 4 (q', -q|Z, I", V) + \/ 4 (-q - q', -q'|Z, I', I") + 
V, (q', q + q'|Z", Z, J') + V 4 (q, q + q'|Z", J', Z)] . (B13) 



APPENDIX C: FOUR-MAGNON AMPLITUDE FACTORS 



The expressions for the amplitude factor Aj (i = 1,2,3) in Eq. (39), and hence for the 
vertices in the diagrams in Fig. 6, are given by 



IX 



A 1 (kk'q|Z 1 Z 2 Z 3 Z 4 )= £ {4fi(k) [SJUCkO^Ufa - k')S n ,, 3 (q - k)5„^(k) + 2^(k') 

nn' 

Tl h (k> - q)T n ,, 3 (k - q)S <u (k) + T n * J-k')T;jk - q)T„,, 3 (k' - q) x 
T n , M (-k) + 2T n * J-k')^ 2 (q - k')S n , h (q - k)T n , u (-k)] + <„,(k') x 
T* 5ii (-k')T^ 2 (k - q)T„, i3 (k' - q)T Bil4 (-k) + 2T n *, J-k')x 
^ 2 (q - k')^ )i3 (q - k)T„, 4 (-k) + 5;, Jk')S n ,Jq - k') x 
S„,, 3 (q - k)S n , k (k) + 2S' n „ f Jk')T;jk' - q)T n ,, 3 (k - q^M - 
2<i(k - k') [5;, Jk')^ 2 (q - k')S n , l3 (q - k)S n ,, k (k)+ 
T:,, h (-k')Tl l2 (k - q)T n , Js (k' - q)T n , M (-k) + S* n , M {k!)T* nM {k! - q) x 
T B , l3 (k - q)5 nV4 (k) + T* , a-k')^ 2 (q - k')^ 3 (q - k)T n , u (-k) 
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2< ( J(q) IS^CkW.Jk' - q)^ 3 (q - k)T n , u (-k) + T n *J-k')x 



SU(q - JOWk - q)5 B ^(k)| + <i,(k) lr;j-k')r;jk' - q 



T n ,< 3 (k - q)5 nV4 (k) + 25;jk')T n *Jk' - q)5 nA (q - k)5 B ^(k) 



+ 



<k"k') |T: Vl (-k')5;jq - k')^ 3 (q - k)S n , u (k) + 2T n *, J-k') 



x 



^ 2 (k' - q)5 nA (q - k)T n , u (-k)j + C^,(k') [^, Il (k')^ Ja (k / - q) x 
T n ,< 3 (k - q)T n , u (-k) + 2^;, h (k')^ 2 (q - k')T n , Js (k - q)S nM (k)] + 

0" k ) Kh(k')^ 2 (q " k')^ 3 (q - k)T„, u (-k)+ 

2^,Jk')T; i2 (k' - q)T rM3 (k - q)T nV4 (-k)] } , (CI) 

A 2 (kk'q|/ 1 / 2 / 3 /4)= E «Uk) [T n * J-q)T n ,, 2 (k + k' - q)T n , Js (-k') W*) + 2^(q)x 

n,n' 

T n ,, 2 (k + k' - q)S nM (k')S n , M (k) + S* njll (q)S nM (q - k - k')S n ,, 3 (k') x 
T„, u (-k) + 2T; i4 (-q)T n ,, 2 (k + k' - q)^, i3 (k')T nV4 (-k)] + 4£,(q) X 
T n *, a-q)T n ,, 2 (k + k' - q)^ )i3 (k')T n , i4 (-k) + ^, JqjT^Jk + k' - q)x 
5 n ,, 3 (k')^, k (k)] - 2< ( J(k - q) [S* n!j Jq)T n ,, 2 (k + k' - q)5 B>l ,(k')+ 
S <k (k) + T*, h (—q)T n j 2 (k + k' - q)^, i3 (k')T„, i4 (-k); - 
2< ( J(k + k') [S* nt aq)^ 2 (q - k - k')T nV3 (-k')S nV4 (k)+ 
T;a-q)T ni , 2 (k + k'-q)T n , i3 (-k')^ 4 (k)] +5jJ,(k) [^(q) x 
S n ,, 2 (q - k - k')S n ,, 3 (k')S <i4 (k) + 2T n *J-q)S n ,, 2 (q - k - k') 
T„, 3 (-k')^,u(k)] + C%l(-k) [r; h (-q)^ 2 (k + k' - q) 
T n , h (-k')T n , M (-k) + 2^ 1 (q)T n ,, 2 (k + k' - q)^, i3 (k') x 
T„, u (-k)] + <i,(-q) [T r t,J-q)5 n , 2 (q - k - k')5 nA (k') 
T n , u (-k)] + C#>,(q) [^(q)^ (k + k' - g) T n , h (-k')S nM (k)} } , (C2) 



Ix 



X 



A 3 (kk'q|/ 1 / 2 / 3 /4)= E «U-k) pO, a-k)^ 2 (k')S;, 3 (q - k - k')S n , k (q)+ 

n,n' 

2T:, h (-k)^ 2 (k')T; i3 (k + k' - q)T rM4 (-q) + S* ^(k) x 
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^ 2 (-k')T; i3 (k + k'-q)T„, 4 (-q) +25;, /i (k)^ 2 (k') x 

T;jk + k' - q)s nM ( q )] + <Uq) K h (Vs: th (k')x 

T^k + k' - q)S nV4 (q) + T;j-k)T; i2 (-k')^, 3 (q - k - k') x 
T n , u (-q)] - 24$ (k - q) [^,ak)^ 2 (k')T^ 3 (k + k' - q)x 
5n',u(q) + a-k)T n ^(-k')^ 3 (q - k - k')T„, u (-q); - 
2< ( J(k + k') [T n *, Zl (-k)^, , 2 (k')^ 3 (q - k - k')S„,, 4 (q)+ 
T^^-k)^ 2 (k')T; i3 (k + k'-q)T„, 4 (-q)] + <>,(-k) X 
T* ;il (-k)r; z2 (-k')T^ 3 (k + k' - q)T n , u (-q) + 2T* )h (-k)x 
^ 2 (k')T^ 3 (k + k'-q)^ 4 (q)] + C#>,(k) k Vl (k)S;, 2 (k') x 



x 



s;, 3 ( q - k - k')s nM ( q ) + 25;, ii (k)^ 2 (k')r;, 3 (k + k' - q ) 

T n , u (-q)] + sW„(q) [^(k)^ J-k')T n ^(k + k' - q) Wo/ . 
+C£i'(-q) [^J-k)^ 2 (k')S; i3 (q - k - k')T n , u (-q)] } . (C3) 

The terms A 4 and A 5 are given by similar expressions. They will not be quoted here, since 
they do not enter into the expressions for the SW energy shift and damping. 

We can now define the weighting factors 6; (i = a, b, c, d) in Eqs. (40) and (41) as 



a = Axfa, q, q + q'|Z, I', I', I) + A^q, q', q + q'|Z', Z, Z) 
+ Axfa', q, q + q'|Z, J', Z, /') + A^q', q', q + q', |Z', Z, Z, J'), 



(C4) 



6 = [Ax(q, q", c£\l"l'"l'l) + A^q, q' - q", q'|Z'TO)+ 

Ai(q' - q, q", c(\l'Tll') + A^q' - q, q' - q", q'|Z'"Z"ZZ')] x 
[AxCq", q, q'|ZZ'Z"'Z") + Mq", q' - q, q'|Z'ZZ"'Z")+ 

Ai(q' - q", q, q'lZZ'ZV") + Ax(q' - q", q' - q, q'|ZOT)] , (C5) 

9 C = [A 2 (q - q', q", q|ZZ'"ZY) + A 2 (q - q', q' - q", q|ZZ"Z'V) + 
A 2 (q", q - q', q|ZZ'"Z'Z") + A 2 (q", q' - q", q|ZZ'Z w Z") + 
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A 2 (q' - q", q", q\U'l"l'") + A 2 (q' - q", q - q', q\ll"lT)\ x 
[A 3 (q - q', q", q\lTl m l) + A 3 (q - q', q' - q", q\lTl"l)+ 
A 3 (q", q - q', q\l"lTl) + A 3 (q", q' - q", d\lT'l'l) + 
A 3 (q - q", q", q|H'77) + A 3 (q' - q',' q - q', q|Z"77'7)] , 

[A 2 (q, q", q' + q?\l m lTl) + A 2 (q, q' - q, q' + q"|/"7'77) + 
A 2 (q", q, q' + q"|/"77Z") + A 2 (q", q' - q, q' + q"|Z'W) + 
A 2 (q' - q, q, q' + c['\l"Tll') + A 2 (q' - q, q", q' + c['\l"'ll"l')] 
[A 3 (q, q", q' + q!'\ll"lT') + A 3 (q, q' - q, q' + <f\UTl m ) + 
A 3 (q", q, q' + cf\l"llT) + A 3 (q", q - q, q' + q ff |W) + 
A 3 (q' - q, q, q' + 4'\l'll'T) + A 3 (q' - q, q", q' + q"|/7'7/ w )] 
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FIGURES 

FIG. 1. The effects of the next-nearest neighbor exchange interactions on the SW spectrum in 
a 16-layer EuO film: (a) SW frequency against J2/J1 for q x = 0.001 (solid curve) and q x = 0.04 
(dashed curve); (b) SW frequency against q x for J2/J1 = (solid curve), J2/J1 = 0.2 (dashed 
curve), and J2/J1 = 0.3 (dot-dashed curve). 

FIG. 2. The linear SW dispertion relation for a 16-layer GdCi3 film with next-nearest neighbor 
exchange J 2 = 0.25Ji solid lines and J 2 = dashed lines. 

FIG. 3. Diagrammatic representation of (a) the non-interacting Green function and (b) the 
possible types of proper self-energy contributions of the system. 

FIG. 4. Diagrammatic representation of the Dyson series for the renormalized Green function 
in terms of the self-energy parts. 

FIG. 5. Second order diagrams for the self-energy contributions due to three-magnon processes. 

FIG. 6. Self-energy diagrams due to four-magnon processes in (a) first order and (b) second 
order. 

FIG. 7. Three-magnon damping of the lowest SW branch (1) vs reduced wave vector q x a/ir 
for EuO films: (a) Total damping (solid curve) and contributions from some of individual (I', I") 
processes: (1,1) dot-dashed curve, (1,2)+(2,1) dashed curve, and (2,2) dotted curve; (b) Total 
damping for the cases of N = 8 (solid curve), N = 16 (dotted line) for the lowest branch 1, and for 
N = 8 and branch 2 (dashed curve). 

FIG. 8. Three-magnon damping for a GdC^ film with 8 layers for SW branch 1. 

FIG. 9. As in Fig. 8, but for a Fe film. 
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FIG. 10. Temperature dependence of the three-magnon damping for a 8-layer EuO film with 
a fixed value of the wave vector (q x /ir = 0.5), showing (a) splitting processes (solid curve), b) 
confluence processes (dashed curve) . 

FIG. 11. SW energy shift vs reduced wave vector (fca/7r for EuO films with N = 8 (solid curve) 
and N = 16 (dotted curve) for the lowest branch 1, and N = 8 for branch 2 (dashed curve). 

FIG. 12. The four-magnon SW energy shift of the lowest branch vs reduced wave vector q x a/ir 
for an EuO film with 8 layers. 

FIG. 13. Temperature dependence of the SW energy shift compared for three-magnon (solid 
curve) and four-magnon (dashed curve) processes. 
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Fig. 5. 
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Fig. 6b. 
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Fig. 7a. 
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Fig. 7b. 



42 



0.00030 



0.00025 



N 

X 

o 



0.00020 



0.00015 







0.2 



0.4 



0.6 



0.8 



q x a/7c 



Fig. 8. 
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Fig. 9. 
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